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Vacuum solutions of Einstein's field equations have played an important role in our 
understanding of curvature effects and relativistic behavior even though it is obvious that 
real stars do not sit in vacuum but have particle and radiation atmospheres. Quantum effects 
allow atmospheres to be added to classical vacuums; for example black holes are associated 
with atmospheres of Hawking radiation [jlj , . In addition to their intrinsic value as exact 
solutions, vacuum solutions in general relativity are approximate string theory solutions for 
curvature small compared to the Planck scale ||. The intense level of activity in string 
theory has lead to the idea that many of the classic vacuum scenarios, such as the static 
Schwarzschild point mass/black hole, may have atmospheres composed of a fluid or field of 
strings |4[]. One of the well known examples of radiation atmospheres is the Vaidya metric 
||, generated from the vacuum Schwarzschild solution by allowing the Schwarzschild mass 
to be a function of retarded time. The resulting stress-energy content describes outgoing 
short-wavelength photons. 

In this letter we extend the Vaidya metric by allowing the mass to be a function of both 
retarded time and distance along the outgoing null geodesies. The effect of this extension 
is to create two fluids outside the star, the original null fluid and a new fluid composed of 
strings. Given the recent links |J, [|7| between black holes and string theories, this result 
is of interest by itself. It is additionally interesting since our new analytic solutions for the 
mass function allow the metric to be written as a superposition of a string fluid and vacuum 
Schwarzschild. We have thereby extended Xanthopolous superposition ||. 

The string fluid tension and density depend on spatial derivatives of the mass function. 
Assuming a specific model for propagation of the density allows the generation of new 
densities and hence new mass functions. We choose to propagate the density diffusively as 
a particular example of mass transport. 

Our sign conventions are 2A u; [ a p] = A^R 11 ua/3 , R^ v = -R a Mra , and metric signature (+,-,- 
,-). Greek indices range over (0,1,2,3) = (u,r,$,<p). m abbreviates dm/du, m! abbreviates 
dm/dr, and m" represents d 2 m/dr 2 . Overhead carets denote unit vectors. 
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The spacetime metric covering the region exterior to a spherical star is given by 

ds 2 = Adu 2 + 2dudr - r 2 {d$ 2 + sin 2 ^ 2 ) (1) 

where A = 1 — 2m(u,r)/r. Initially m(u,r) = mo provides the vacuum Schwarzschild 
solution in the region r > 2m . At later times m(u,r) admits a two-fluid description of 
diffusing matter and outward flowing short-wavelength photons (sometimes called a "null 
fluid"). Metric (HD is spherically symmetric and given in retarded time coordinate u. With 
the use of a Newman-Penrose null tetrad the Einstein tensor is computed from ([!]) and given 
by 

G^ = -2$n(/^n Jy + n^ v + m^m u + m^m v ) (2) 
-2$22y„ - 6A^. 

Here the null tetrad components of the Ricci tensor are 

$ n = -ij-m" - 2m')/(4r 2 ), (3a) 
$ 22 = -m/r 2 , (3b) 
A = R/24 = (rm" + 2m) /(12r 2 ). (3c) 

The metric is Petrov type D with and principal null vectors, geodesic, and 

l^dx^ = du, (4a) 
n^dx^ = (A/2)du + dr, (4b) 
m^dx^ = —(r/y/2)(d& + i siw&dcp). (4c) 

In order to clearly see the two-fluid description we introduce a timelike unit velocity vector 
and three unit spacelike vectors f M , (p^ such that 

The unit vectors are defined by 
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VpdxP = A^ 2 du + A~ 1/2 dr, (5a) 

r^dx? = A~ 1/2 dr, (5b) 

ti^da? = rd$, (5c) 

if^dx^ = rsmddip. (5d) 

The Einstein tensor (§) can be written as a two-fluid system: 

G„ v = {2fn/r 2 )l^l u - {2m' /r 2 ) {v„v v - r^r u ) (6) 
+ {m"/r){$^ u + fl u ). 

Spherical symmetry allows the function m(u, r) to be identified as the mass within two- 
surfaces of constant u and r, and invariantly defined from the sectional curvature of those 
surfaces: 

- 2m/r 3 = R aPlx J a ^S^. (7) 
The string bivector is defined by 

EP " = £ fs^' (B.C) = (0.1) O r(2,3). 

Spherical symmetry requires that the string bivector have a world-sheet in either the (u, r) 
or (■&, ip) plane. We require that the world-sheets be timelike, i.e. 7 := ^E^E^ < 0, and 
so only the E ur component is non-zero. Here 7 = — 1. It is useful to write YJ* ,V in terms of 
unit vectors 

IT = r»v v - fr", (8) 

and so E^E^ = v fl v u — f fl f u . We follow Letelier ||, [H| and write a string energy- 
momentum tensor by analogy with one for a perfect fluid. The string energy-momentum is 
given by 



Turing = p{ _ l) l/2±«± oiv _ p±Ht 
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where H% = &> v - t» a t av , E^ v t vfi = 0. We have written := (--f)' 1 / 2 ^ , so that 
is invariant under reparameterizations of the world-sheets [g]. Einstein's field equations, 
G /tl/ = —8iiT^ v , allow the matter portion of © to be identified as a string fluid: 

Tiiv = Wnlv + + Vrr il r v + p^tf/A + 0p(p v ). (9) 

Thus 

47r^ = — rh/r 2 , (10a) 

471-p = -47rp r = m'/r 2 , (10b) 

87rp_i_ = —mf'/r. (10°) 

The equation of motion = is identically satisfied for the components of 7) t „ given in 
Eqs.(H). 

As an example of mass transport we assume the strings diffuse and that string diffusion 
is like point particle diffusion where the number density diffuses from higher numbers to 
lower according to 

d u n = V V 2 n. (11) 

V 2 = r~ 2 (d I dr)r 2 (d / Or) , and D is the positive coefficient of self-diffusion. Classical trans- 
port theory derives the diffusion equation by starting with Fick's law 

J (n) = -VVn (12) 

where V is a purely spatial gradient. Then 4-current conservation Jjty^ = 0, where 

= nd u — V(dn/dr)d r , 

yields the diffusion equation ([TT]) . We label the 4-current J( n ) to indicate n diffusion but we 
could have also written Ji p \ since the string number density and string fluid density must 
be related by p = M s n where M s is the constant mass of the string species. M s must be 



a multiple of the Planck mass since it is only over Planck length scales that point particles 
resolve into strings. 

By rewriting the components ( |10a| ) and ( |10b| ) as fa = —47rr 2 ifj and ml = 4irr 2 p, we 



can write the integrability condition for m as 

p + r- 2 d r (r 2 ?p) = 0. (14) 
If we compare the diffusion equation (|Tl|) (n replaced by p) 

p = V r~ 2 d /dr(r 2 dp/dr) (15) 

with p in Eq.(|l~4D we obtain 

m = AnV r 2 dpj dr. (16) 

Thus solving the diffusion equation for p and then integrating those solutions to obtain m 
provides exact Einstein solutions for diffusing string fluids. 

There are many analytic solutions of ( |ToD and three of them are 

p = p + k 1 /r, (17a) 
p = po + k 3 u- 3/2 exp[-r 2 /(4Vu)}, (17b) 
P = Po + (h/r){l + ( v / ^/2)erf[r(4P M )- 1 / 2 ]}. (17c) 

Upon integrating m' = 47rr 2 p and rh = 4tcV r 2 dp/dr we obtain the following masses, 
listed consecutively, for the densities above: 

m(u, r) = m + (4?r/3)r 3 po + 2nk 1 (r 2 - 2Vu). (18a) 
m(u, r) = m + (47r/3)r 3 p + 167rA;3D 3/2 [-r/exp(-r/ 2 ) + (v^F/2)erf(r/)]. (18b) 
m(u, r) = m + (47r/3)r 3 p (18c) 



+2nr 2 h{(l - 2^2) + [(v^/2)(l - ^)<*Kv) + (2r / )- 1 exp(-r / 2 )]}. 



where r\ := r (4Pm) 



'1/2 
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It is clear that metric ([D can be written in Kerr-Schild form as r]^ v — [2m(«,r)/r]y„. 
When m(u,r) = m for the Schwarzschild solution, (2mo/r)l^l u solves the vacuum field 
equations linearized about flat space. This was Xanthopoulos' original superposition. His 
generalization || has g^ v = + if (x^Z^,, with if(x M )Z M / v required to solve the vacuum 
field equations linearized about g^. Referring to the mass solutions above, we can write 
1 — 2m(u,r)/r = 1 — 2mo/r — 2m/r. Metric (H) then has the form rj^ — (2mo/r)l^l v — 
(2fh/r)l^l v which is clearly 

9r = C h - (2m/r)y, (19) 

For the Vaidya metric, with m(-u), the field equations linearized about g^ h yield 
GjJ) = (dm/ ' du)r~ 2 l^l v |TT| which is not a vacuum solution and so Gfy computed about 
Schwarzschild with m(u, r) is a fortiori not zero. Since g^ u is an exact solution for a string 
fluid and g^ h is an exact vacuum solution, we have extended Xanthopolous' generalization. 

E.N. Glass was partially supported by an NSERC of Canada grant. Computations were 
verified using MapleV.4 (Waterloo Maple Software, Waterloo, Ontario) and GRTensorll rel 
1.59 (P. Musgrave, D.Pollney, and K. Lake, Queens University, Kingston, Ontario). 
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